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The Absolute Abundance of Iron in the Solar Corona 
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ABSTRACT 

We present a measurement of the abundance of Fe relative to H in the solar 
corona using a technique which differs from previous spectroscopic and solar wind 
measurements. Our method combines EUV line data from the CDS spectrometer 
on SOHO with thermal bremsstrahlung radio data from the VLA. The coronal 
Fe abundance is derived by equating the thermal bremsstrahlung radio emission 
calculated from the EUV Fe line data to that observed with the VLA, treating 
the Fe/H abundance as the sole unknown. We apply this technique to a compact 
cool active region and find Fe/H = 1.56 x 10~^, or about 4 times its value in the 
solar photosphere. Uncertainties in the CDS radiometric calibration, the VLA 
intensity measurements, the atomic parameters, and the assumptions made in 
the spectral analysis yield net uncertainties ~ 20%. This result implies that low 
first ionization potential elements such as Fe are enhanced in the solar corona 
relative to photospheric values. 

Subject headings: atomic data — Sun: abundances, atmosphere, corona, radio 
radiation. X-rays 



1. INTRODUCTION 



It seems remarkable that the absolute abundances (here taken to mean abundances 
relative to hydrogen) of trace elements in the solar corona are still a matter for debate when 
many measurements for much more distant stars exist. The reason for this situation is easy 
to understand. Classic spectroscopic techniques for measuring relative abundances require 
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that we compare the amount of radiation emitted in spectral hnes whose atomic physics we 
understand and which arise from the two different elements in question. In the solar corona 
H and He are fully ionized, and hence produce no lines. In stars with hotter coronae, thermal 
bremsstrahlung produces a continuum emission which can be used instead of a spectral line 
as a measure of the amount of H present, but the Sun's quiescent corona is so cool that the 
continuum present at X-ray or EUV wavelengths is relatively weak and difficult to measure. 

It has been known for some time that there is a difference in the relative abundances 
of elements in the photosphere and the corona organized according to the first ionization 
potential (FIP) of the element: high-FIP elements such as O are underabundant in the 
corona relative to low-FIP elements such as Fe by a factor of 4 (Meyer 1985). However, 
it has not been clear whether low-FIP elements are enhanced in the corona in an absolute 
sense (i.e., as a fraction of the total ion number density, dominated by H and He) or the 
high-FIP elements are depleted relative to H. Two main techniques have been used to infer 
absolute abundances in the solar corona (Fludra et al. 1999). One can observe solar flares, 
which produce plasma hot enough for the continuum to be measured, and compare the 
continuum with a suitable set of lines; or one can measure particles directly in the solar 
wind, and assume that the abundances there truly reflect the abundances in the solar corona. 
Unfortunately, these two techniques have resulted in conflicting measurements for, e.g., the 
absolute abundance of iron. The flare measurements suggest that it is close to the absolute 
abundance in the photosphere (Veck & Parkinson 1981) or perhaps factors of 1.5 - 2 larger 
(Fludra & Schmelz 1999), while the solar wind particle measurements suggest that it is 4 
times the photospheric ratio (Reames 1992). 

In this paper we use a third technique to measure the absolute abundance Npe/Nn 
which has several desirable features. We use EUV observations of lines of Fe in order to 
determine the amount of Fe present in the corona, but to measure the amount of H present 
we use radio continuum measurements of thermal bremsstrahlung from coronal plasma. This 
measurement therefore pertains to the quiescent solar corona, including non-flaring weak 
active regions, rather than flare or solar wind plasma. Like classic spectroscopic techniques, 
this technique should be independent of the fllling factor of plasma in the volume studied 
since both types of emission have the same dependence on number density. We describe 
the measurement technique and the observations used for this measurement in the following 
section. We then discuss the implications of our measurement. 



-3- 



2. THE MEASUREMENT TECHNIQUE 

The method proposed here to measure the abundances of elements relative to H rests on 
the fact that the EUV line emission and radio bremsstrahlung both depend on the product 
of the electron number density and the number density of an element. In the case of radio 
bremsstrahlung the element is effectively hydrogen; in the case of an EUV line, it is the 
element responsible for the hne. Therefore comparison of the two yields the abundance of 
the latter element relative to hydrogen. 

Specifically, since most EUV lines produced in the solar corona are coUisionally ex- 
cited and optically thin, the line intensity at a photon energy E corresponding to a tran- 
sition in ionization state Z of element X is the integral over temperature of the product 
Qx,z,eDEM{T)Nx/Nu (Fludra et al. 1999), where Qx,z,e is the "contribution function": 
it includes the intrinsic spontaneous emission rate (the Einstein "A" coefficient) as well as 
the ratio Nx,z/Nx of the number density of the particular charge state Z of the element X 
responsible for the line to the total number density of the element. This ratio is usually de- 
termined by an ionization balance calculation appropriate to the temperature of the plasma, 
and Q is a strong function of temperature for this reason. However, for most lines Q is 
effectively independent of Nx and Ng. The (hydrogen column) differential emission measure 
DEM(T) is the integral along the hne of sight of d{Ne N^^j/dT. 

Bremsstrahlung is one of the two main mechanisms for radio emission in the solar corona 
and also has a significant dependence on temperature: in the optically thin limit, the fiux 
at a frequency / is 

5=9.78x10-3^(1 + 4^) / / T-^-'DEM(T)G(T) dTdn (1) 

Nh J J 

in ergs cm'^ s"^ Hz"^ (Bulk 1985), where G{T) = 24.5 - ln(T//) is the Gaunt factor, ks 
is the Boltzmann constant, c the speed of fight and the second integral is over the solid 
angle subtended by the source. We are assuming that in a plasma at the temperature of 
the solar corona the ion component is completely dominated by protons and fully-ionized 
hefium, and we adopt Nue/Nn — 0.1. It is clear from (1) that we must know the DEM in 
order to calculate the thermal bremsstrahlung radio flux accurately. Note that as long as 
Nx/Nu does not vary from place to place, (1) has exactly the same dependence on density 
as the EUV line emissivity and thus filling factor has exactly the same effect on the EUV 
line and radio continuum fiuxes. 

We determine the abundance A^Fe/-^H with the following steps: (i) We measure the 
radio fiux from a source in the solar corona. This source must be optically thin and its 
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radio flux must be due only to bremsstrahlung. The latter requirement can be tested by 
observing the radio spectrum of the source, which should be essentially flat for optically thin 
bremsstrahlung. (ii) We observe the same region with the Coronal Diagnostic Spectrometer 
(CDS) on SOHO in order to measure the intensities of lines from a range of charge states of 
Fe sufficient to give a reliable determination of the DEM. (iii) We use the CHIANTI atomic 
data package (Dere et al. 1997; Landi et al. 1999) to fit the Fe hue fluxes to a DEM. (iv) We 
assume that the shape of the DEM for electrons and protons is the same as that for Fe, and 
integrate over temperature in (1) to predict the optically-thin bremsstrahlung radio flux we 
expect to observe, for a given absolute abundance of Fe. Comparison with the radio flux 
measured by the VLA yields Npe/N-n. 

Coordinated observations with CDS and the VLA^ were carried out on 1997 July 27, 
Aug. 3 and Nov. 11. We will present an analysis of all the data elsewhere; here we discuss 
the best target observed, a small active region AR 8105 found at about S35W15 on 1997 
Nov. 11 (with the VLA in "D" configuration). This was the best target in the sense that it 
was compact: the VLA is an interferometer and cannot measure the flux of large sources, so 
radio flux measurements are more reliable for compact sources. The region was one of four 
observed on this day with the VLA at 1.4, 4.8 and 8.4 GHz. CDS also observed the same 
4 regions alternately using the sequence FEJONS which reads nine 21-pixel-wide windows 
from the NISI detector of CDS containing 14 identified lines including those of Fe X, XI, XII, 
XIII, XIV, XV, XVI and XVII. One run of the sequence took 32 minutes, scanning with the 
4" sht across a 240" wide field of view. Five sequences were acquired on AR 8105 centered at 
about 13:14, 16:28, 18:42, 20:56 and 23:10 UT, respectively. The VLA observations covered 
the period from 14:30 UT to 23:10 UT, cycling between the different frequencies rapidly in 
order to achieve optimal uv coverage. Images of AR 8105 acquired by CDS, the VLA and 
EIT are shown in Figure 1. 



3. THE ABUNDANCE OF IRON 

It is essential for this measurement that the VLA see the same plasma that is seen 
in the EUV, and Fig. 1 shows that this is indeed the case: the EUV and radio images 
match each other extremely well, indicating that gyroresonance opacity docs not contribute 
signiflcantly to the radio emission and we conclude that it is due to bremsstrahlung. The 
higher-resolution EIT images shows that two east-west loops dominated AR 8105. 
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The VLA observations resulted in images with resolutions of 30" at 1.4 GHz, 13" at 4.8 
GHz and 7" at 8.4 GHz. Using the 4.8 and 8.4 GHz images we identified an area of dimension 
69" X 97" around AR 8105 which includes essentially all the fiux from the region. Radio 
fluxes were derived by summing over this area at 4.8 and 8.4 GHz, while due to the poorer 
resolution at 1.4 GHz we instead fitted a Gaussian model to the feature associated with AR 
8105 and determined its fiux. At each frequency we have independent measurements in two 
separate 50 MHz bands and in the two circular polarizations, all of which are consistent 
with one another. The resulting fluxes are 0.25 ± 0.01 sfu at 1.4 GHz, 0.207 ± 0.004 sfu 
at 4.8 GHz and 0.197 ± 0.006 sfu at 8.4 GHz, which are consistent with the /~°'^ spectrum 
expected for optically thin bremsstrahlung due to the Gaunt factor in (1). The quoted 
uncertainties are the map rms per beam times the square root of the number of beams over 
which integration was carried out, and are much larger than the formal uncertainties in the 
fluxes. Polarizations were also low, consistent with the properties of bremsstrahlung. 

The GDS data were calibrated using the standard reduction package in SOLARSOFT 
(results presented here use the version current in 1999 December). Line proflles were fltted 
and fluxes extracted, with uncertainties calculated according to Thompson (1998). We used 
both a recent radiometric CDS calibration derived from comparison with the SERTS rocket 
flight in 1997 (Thomas et al. 2000) and the "version 2" CDS calibration which was in the 
CDS software as of 1999 December. One important step is that background subtraction was 
carried out on the CDS data using a hnear fit to pixels to the north and south of the region 
studied, well away from AR 8105. This step is important because of the fact that the VLA 
does not measure large-scale flux and hence is insensitive to the quiet-Sun contribution. In 
effect we assume that a smooth quiet-Sun contribution has been subtracted from the radio 
fluxes, and hence it should be subtracted from the CDS data as well. This predominantly 
affects the cooler lines, amounting to up to 20% of the photons in the original spectrum, but 
being fairly uniform across the CDS fleld of view. 

The following lines (with typical line fluxes averaged over the region shown in paren- 
theses, in ergs cm~^ s~"^ st~^) were then used to calculate the DEM: Fe X at 345.7 A (55 ± 
3), Fc XI at 352.7 A (140 ± 4), Fc Xll at 352.1 A (145 ± 4) and 364.4 A (230 ± 4), Fe XllI 
at 348.3 A (95 ± 4), Fe XIV at 334.2 A (260 ± 4), Fe XVI at 360.7 A (380 ± 4) and Fe 
XVII (undetected here) at 350.5 A. These fines were chosen to be as density-insensitive as 
possible while still covering the temperature range expected for the non-flaring solar corona, 
log T = 5.8 up to 6.8. The Fe XIII and Fe XII fines do have significant sensitivity to density 
in the likely range of coronal densities, 10^ - 10^° cm~^. In practice, AR 8105 was quiet and 
there is no significant plasma outside the temperature range 5.9 - 6.4. The fitting procedure 
uses bins of width 0.1 in log T; in effect, we find that only 4 bins contain significant emission 
measure and the adjacent two bins are also significant, so we are fitting to 6 parameters using 
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8 lines from 7 different cfiarge states. The Fe XV line at 327.1 A was not used to determine 

the DEM because it was always found to be ~ 30% lower than the models predicted; the 
atomic data for it are clearly not consistent with the strong Fe XIV and XVI lines observed. 
Additional details will be given in a companion paper. 

We adopted the ionization equilibrium calculations of Arnaud & Raymond (1992) and 
assume a constant pressure nT of 1 x 10^^ cm~^ K (e.g., = 5 x 10^ cm~^ at T = 2 x 
10® K) in determining the DEM. The line fluxes used for the fit were generally reproduced 
by the resulting model DEM to better than 10% (SERTS) or 5% ("version 2"). Figure 
2 shows the geometric mean DEM over the 5 CDS sequences with error bars showing the 
standard deviation in each temperature bin. Essentially no plasma is present below log T 
= 6.0 (controlled by the ratio of Fe X to Fe XI) while the amount of plasma above 2 x 10® 
K is small (controlled by the ratio of Fe XVI to Fe XIV). These properties are consistent 
with those of DEMs derived for other regions of similarly low activity (Brosius et al. 1996). 
Integrating over each of the 5 DEMs resulted in a predicted radio flux of 0.829 ± 0.023 sfu 
in the 6700 arcsec^ area at 4.8 GHz if we use the CDS calibration based on the SERTS- 
97 comparison and adopt Npe/Nn = 3.9 x 10~^ as in Meyer (1985), which is close to the 
photospheric value. The observed radio flux of 0.207 sfu is a factor of 4.0 smaller than 
this. Consequently we derive a value of 1.56 x 10^'^ for ApeZ-^H, or log^Q /Ipe = 8.19 on 
the commonly used scale where logj^g^H ~ 12.00. Using the "version 2" CDS calibration 
we predict 0.716 ± 0.013 sfu and Npe/Nu = 1.35 x 10"^ (logio ^Fe = 8-13)- We adopt 
the more recent SERTS-97 cahbration for our final result. The uncertainty of 3% for the 
predicted radio fiux is the standard error in the mean of the results of the five distinct CDS 
observations of AR 8105 and thus incorporates the effects of both temporal variability over 
the 12 hours of the experiment as well as the stochastic nature of the process of fitting the 
line fiuxcs to a DEM. We adopt an uncertainty of 15% for the CDS radiometric calibration 
(Thomas et al. 2000). 

It turns out that the measured value for Npc/Nu is reasonably robust to several of the 
assumptions made here, mostly because it depends on an integration over the DEM and is 
thus not very sensitive to details of the DEM. One assumption is that the result quoted uses 
the version of CHIANTI (0.9) used by the CDS team at NASA/GSFC as of 1999 December, 
which is not the most current version (2.0). For the Fe fines used here the only differences 
between the two versions appear to be in Fe X and Fe XI, and the newer version results in 
significantly poorer fits for the DEM because of these two lines: the resulting DEMs have 
less emission measure at low temperatures and more at high temperatures. Notwithstanding 
this fact, the result for the abundance differs by less than 3% with the newer version of 
CHIANTI. If we use the ionization equilibrium calculations of Arnaud & Rothenflug (1985), 
the DEMs show a pronounced peak at logT = 6.3 and the resulting abundance is about 
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10% larger. A change in the assumed pressure by a factor of 2 affects the results by a 
similar amount: it would be 8% smaller at 5.0 x 10^^ cm~^ K and 10% larger at 2.0 x 
10^^ cm~^ K. We attempt to minimize the uncertainty due to variation in pressure across 
the active region by choosing a value representative of the source region for the bulk of the 
EUV line emission. We note also that the background subtraction carried out on the CDS 
data introduces some uncertainty in the final results, as does solar calibration at the VLA, 
which uses dedicated noise sources for solar calibration which need to be measured routinely 
to maintain calibration accuracy. However, the calibration at each of the three frequencies 
observed is independent of the others, so we take the fact that the measured spectrum is 
flat to be an indication that solar calibration was not in error by more than the 20% levels 
which all the other factors may introduce. 



4. CONCLUSIONS 

Our measurement of an absolute abundance for Fe which is 4 times photospheric is 
clearly consistent with the interpretation that low-FIP elements such as Fe are enhanced 
relative to their photospheric abundances. Of the commonly-used solar abundance tables, 
Feldman's (1992) is closest to our measurement with log^g ^Fe = 8.10; his value is based more 
on the argument that H should be treated as a high-FlP element rather than any specific 
measurement of A^pc/A^H- Other commonly used values for log^Q Apc arc the photospheric 
results 7.60 (Allen 1973), 7.51 (Anders & Grevesse 1989) and 7.67 (Grcvcsse & Anders 1991), 
the coronal results 7.59 (Meyer 1985), 7.83 (Fludra & Schmelz 1999) and 8.50 (Waljeski 
et al. 1994), and the solar wind measurements 7.93 (Reames 1995) and 8.51 (Reames 1999). 
We note in passing that spectroscopic measurements of the Fe abundance in active stars 
with coronae much hotter than the Sun's have tended to suggest abundances smaller than 
photospheric (Mewe et al. 1996; Drake 1996), in contrast to the solar behaviour we find here. 

We have reported here a single measurement of the Fe abundance. Spectroscopic studies 
of relative abundances suggest that they may vary considerably from flare to flare and from 
one region to the next (Sylwester et al. 1984; Saba 1995; Sylwester et al. 1998), so it is 
clearly important to conflrm our measurement. We have additional data sets which will be 
published in a longer paper elsewhere, but with larger uncertainties expected in the radio 
fluxes, and expect to carry out further studies of this type in the future. 
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Fig. 1. — Images of the target region used for this study obtained by CDS (left panel, all 
photons, smoothed with a 7" gaussian), VLA (second panel, 4.8 GHz, beam size 13"), EIT 
Fe XII 195 A (third panel, greyscale, 5'.'2 resolution) and the EIT image smoothed to 10" 

resolution and contoured for direct comparison with the CDS and VLA images. The CDS 
and EIT images have been background-subtracted as described in the text. In the contour 
images, contours are plotted at 5%, 10%, 95% of the maximum for that image. 
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Fig. 2. — The geometric mean of the best-fit differential emission measure distributions for 
the five CDS observations of the region studied. The sohd hne shows the resuhs for the 
SERTS 1997 calibration of CDS, while the dotted line shows the results for the "version 2" 
calibration. Error bars show the standard deviation (in the log) over the five measurements 
in each temperature bin. 



